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The  objective  of  the  present  work  was  to  experimentally 
investigate  the  effect  of  small  sinusoidal  perturbations  on 
the  breakdown  of  circular  water  jets  issuing  from  long  tubes 
aiming  to  examine  how  the  mean  velocity  profile  for  a  fully 
developed  laminar  ,  transitional  and  turbulent  flow  as  well 
as  the  frequency  affects  the  mechanism  of  breakdown  of  the 
jet  f 1 ow . 

Observations  using  a  high  frequency  stroboscope  were 
recorded  using  a  high  speed  photographic  technique  ,  for  the 
most  distinctive  phenomena  observed,  for  three  selected  Rey- 
nolds numbers  corresponding  to  laminar  ,  transitional  and 
turbulent  flow  .  Markedly  different  behavior  of  flow  break- 
down and  droplet  formation  was  observed  for  each  of  these 
flows  as  frequency  waveform  or  phase  angle  was  altered. 
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I .   INTRODUCTION 

An  important  physical  phenomenon   which  has  many  practi- 
cal consequences  in  technology,  ranging  from  fuel  injection 
sprays  to  ink  jet  printers,  is  the  breakdown   of  a  circular 
jet  of  fluid  issuing  into  a  quiescent  environment  . 

This  problem  has  been  extensively  studied  previously 
by  Rayleigh  and  others,  [Refs.  1,  2  and  3],  for  steady  uni- 
form jets  issuing  from  an  orifice  ,  but  only  recently  have 
the  effects  of  jet  velocity  profiles  or  other  initial  condi- 
tions on  the  flow  been  investigated,  [Refs.  4,  5  and  6]. 

Theoretical   investigators  have   typically   assumed  jet 
breakdown  to  be  physically   similar  to  transition  from  lami- 
nar to  turbulent  flow   and  have  examined  the  phenomena  using 
classical  stability  theory. 

Rayleigh  ,  [Ref.  1],  studied  the  inviscid  stability  of 
Poisseuille  flow  and  concluded  that  the  flow  was  stable  to 
infinitesimal  disturbances.  This  work  attracted  the  interest 
of  many  investigators  and  sixty  years  later  Tatsumi  [Ref.  2] 
published  his  work  related  to  the  study  of  the  velocity  dis- 
tribution in  the  axisymmetric  laminar  inlet  flow  through  a 
circular  tube,  under  the  assumption  of  '  almost  similarity  ' 
of  velocity  profiles  .  He  obtained  fairly  good  agreement  of 
theoretical  with   experimental  velocity   distribution  curves 
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in  a  rather  limited   region  near  the   entrance  where  the  ap- 
proximation was  valid  . 

Batchelor  and  Gill,  [Ref.  3],  have  analyzed  the  stabili- 
ty of  steady  axisymmetric  jets  considering  the  typical  small 
disturbance  as  a  Fourier  component  with  sinusoidal  dependen- 
ce on  both  axes  in  cylindrical  coordinates  . 

Crow  and  Champagne,  [Ref.  4]  ,  have   employed   a  loud- 
speaker  to   generate  a  wave  train   on  a  turbulent   jet  and 
found  that  the  phase  velocity  of  the  waves  could  be  descri- 
bed  by  applying  linear  theory  of  temporally  growing  insta- 
bility to  a  'top  hat'  velocity  profile. 

Salwen  ,  Chester  and  Grosch,  [Ref.  5],  have  studied  the 
stability  of  Poisseuille  flow  in  a  pipe  of  circular  cross 
section  with  axisymmetric  disturbances  .  They  formulated  the 
problem  by  expanding  the  perturbation  velocity  and  pressure 
in  a  complete  set  of  orthonormal  functions  which  satisfied 
the  boundary  conditions  . 

In  a  recent  work  by  Anderson  and  Bejan,  [Ref.  6],  the 
authors  presented  a  linear  stability  analysis  of  the  large 
scale  structure  of  a  round  jet  surrounded  by  an  annular  she- 
ar layer  .  Their  study  was  limited  to  the  developing  region 
near  the  jet  nozzle  at  very  large  Reynolds  number  (  Re-*-<*>  )  . 
They  examined  the  radial  dependence  of  the  amplitude  of  gro- 
wing disturbances  in  order  to  illustrate  the  extent  to  which 
the  disturbance  penetrates  into  the  jet  and  its  surroundings 
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and  found  that  the  region  influenced  by  a  disturbance  was 
directly  proportional  to  the  wavelength  of  the  disturbance. 
They  also  pointed  out  that  amplified  disturbances  exhibit  a 
phase  lag  across  the  shear  layer  ,  which  may  account  for  the 
spade-like  structures  evident  in  flow  visualizations  of  tur- 
bulent jets. 

In  the  present  work  we  seek  to  investigate  experimental- 
ly the   effects  of   small   sinusoidal   perturbations  on  the 
breakdown  of  circular  jets   issuing   from  long  tubes  ,  which 
permit  some   control  of  their   mean  velocity  profiles  ,  for 
both  fully  developed  laminar  and  turbulent  flows. 
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II.   EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

A.   DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 

The  experiment  was  carried  out  in  a  vertical  falling 
jet  produced  at  the  outlet  of  a  70.0  inch  long  length  of 
0.1875  inch  inside  diameter  stainless  steel  tube  .  The  inlet 
of  the  tube  was  fed  from  a  phenum  chamber  machined  into  an 
aluminum  block  ,  which  also  contained  a  75  Watt  loudspeaker 
driver  unit  mounted  properly  to  this  block,  used  to  introdu- 
ce the  pressure  perturbations  .  A  schematic  diagram  of  the 
experimental  apparatus  is  illustrated  in  Fig.  1  .  Details  of 
the  aluminum  block  and  adapter  are  shown  in  Fig.  2. 

Pressure  fluctuations  as  well  as  mean  pressure  within 
the  phenum  chamber  was  monitored  by  the  Statham  differential 
pressure  transducer  shown  in  Fig.  1  and  a  mercury  manometer. 
Mean  water  flow  was  controlled  with  a  needle  valve  and  mea- 
sured with  a  0.6  gpm  rotameter  as  shown  in  Fig.  1.  In  order 
to  minimize  effects  of  geometry  at  the  exit  end  of  the  steel 
tube,  it  was  very  carefully  machined  to  ensure  sharp  edges 
and  an  absence  of  burrs. 

The  loudspeaker  driver  was  excited  with  a  Wavetec  audio 
oscillator  capable  of  producing  sine  waves  ,  as  well  as 
triangular,  unit  functions  and  ramp  functions  from  1  Hz  to 
10  MHz,  amplified  by  a  Hewlett  Packard  power  amplifier. 
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Wavetec  H.F  VCG  Generator  ,  Model  142 

Hewlett  Packard  Power  Amplifier  ,  Model  467A 

Tektronix  Dual  Beam  Oscilloscope  ,  Type  551 

High  Frequency  Stroboscope  ,  Model  839 

Cimcon  Digital  Multimeter  DMM  51 

Hewlett  Packard  RMS  Voltmeter  ,  Model  3400A 

Intronics  Power  Supply  +15  VDC 

Calibration  Bridge 

110  VAC  Outlet 

Rotameter  0.6  gpm 

Statham  Pressure  Transducer  PM131  TC 

University  Sound  Driver  Unit  ID-75 

Mercury  Manometer 

Aluminum  Block 

Needle  Valve 

L.E.D  Strobol ight 


Figure  1.   Schematic  Diagram  of  Experimental  Apparatus 
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Excitation  voltage  as  well  as  the  output  of  the  plenum  cham- 
ber pressure  transducer  were  monitored  on  a  two  channel  Tek- 
tronix oscilloscope. 

Visual  observations  of  the  jet  breakdown  process  was 
made  with  the  aid  of  a  high  frequency  LED  stroboscope,  deve- 
loped at  the  Max-Planck  Institute  for  Fluids  Research  which 
produces  a  constant  mean  light  output  from  100  Hz  to  1  MHz. 
This  device  not  only  may  be  set  for  any  phase  angle  but  also 
contains  circuitry  which  produces  automatic  phase  angle 
sweep  at  an  adjustable  rate  which  permits  periodic  phenomena 
to  be  visualized  as  they  develop.  Unfortunately  the  light 
output  was  not  sufficient  for  photographic  recording. 

Photographic  recording  of  the  observed  phenomena  was 
made  with  the  illumination  from  a  General  Radio  Microflash 
unit,  which  altough  brighter  than  the  LED  stroboscope,  lac- 
ked the  sophisticated  triggering  circuitry  to  permit  control 
of  the  phase  angle  .  A  4x5  inch  view  camera  fitted  with  a 
210  mm  objective  was  used  to  photograph  on  Kodak  Tri-X  film 
which  was  force  developed  to  achieve  an  effective  film  speed 
of  ISO  3200  . 

An  overall  view  of  the   experimental  apparatus  concern- 
ing and  the  photographic  equipment  used  is  shown  in  Fig.  3. 


xThe  author  extends  his  gratitude  to  Dr.  G.E.A  Meier 
who  made  the  instrument  available  for  this  work. 
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Figure  3.   Photo  of  Overall  Experiment  Setup 
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B.   EXPERIMENTAL  PROCEDURE 

The  first  part  of  the  experiment  was  devoted  to  calibra- 
tion and  determination  of  the  steady  flow  behavior  of  the 
jet  without  the  introduction  of  external  perturbations. 
As  the  flow  rate  was  increased  a  laminar  flow  column  was  ob- 
served to  develope  which  converged  and  began  to  breakdown 
into  a  turbulent  flow  at  100  to  125  diameters  downstream  of 
of  the  tube  exit  corresponding  to  Reynolds  numbers  of  880  to 
2200  . 

Increasing  the  Reynolds  number  above  2200  resulted  in 
a  turbulent  jet  which  was  sudden  and  easily  identified  with 
the  naked  eye  .  Breakdown  of  the  turbulent  jet  was  observed 
with  the  high  frequency  stroboscope  and  was  observed  to  ran- 
ge from  55  diameters  downstream  at  Reynolds  number  of  2288 
to  25  diameters  as  the  Reynolds  number  was  increased  to  4400 

Based  on  these  preliminary  results  ,  a  subcritical  Rey- 
nolds number  of  1760  was  adopted  for  laminar  flows  and  2640 
for  turbulent  flows  and  the  effects  of  small  perturbations 
on  the  flow  at  these  two  Reynolds  numbers  was  investigated 
over  the  frequency  range  of  0  to  2000  Hz  using  the  high 
frequency  stroboscope  to  observe  the  flow. 

Based  on  these  observations  ,  photographic  records  were 
made  at  selected  frequencies  at  which   significant  characte- 
ristic  phenomena   occured  .  Because   of   the   difficulty  in 
accurately   describing  the   pressure   fluctuations   actually 
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produced  in  the  phenum  chamber,  oscilloscope  traces  of  these 
pressure  variations  as  well  as  the  driving  signal  were  also 
recorded.  These  results  are  presented  and  discussed  in  Chap- 
ter IV.  A  discussion  of  the  effects  of  plenum  pressure  vari- 
ation on  the  flow  at  the  exit  of  a  long  tube  is  presented  in 
Chapter  III  .  An  analysis  of  the  dynamic  response  of  the 
pressure  transducer  and  its  calibration  curves  may  be  found 
in  Appendices  C  and  D. 
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Ill .  THEORETICAL  ANALYSIS 

A.   BREAKDOWN  OF  LAMINAR  PIPE  FLOW 

One  of  the  most  important   problems  in  fluid  mechanics, 
one   which   has  attracted   the  interest  of  investigators  for 
many  years  ,  is  that  of  the  transition  from  laminar   to  tur- 
bulent flow. 

The  theoretical  analysis   of  the  stability   of  laminar 
flow  is  based   upon  the  assumption  that  turbulence  is  produ- 
ced by  the  growth   of  what  are  initially  small  perturbations 
or  disturbances  of  the  basic  ,  laminar  velocity  profile. 

The  theoretical  analyses  have  generally  been  limited 
to  infinitesimal  disturbances  for  which  the  differential 
equations  of  the  disturbance  can  be  linearized. 

The  basic  problem  of  the  small-disturbance  theory  is 
to  determine  whether  a  small  disturbance  is  amplified  or 
damped.  If  any  small  disturbance  is  amplified,  it  is  assumed 
that  a  transition  to  turbulent  flow  will  ultimately  occur. 
However  ,  the  small  disturbance  theory  cannot  reveal  any  of 
the  details  of  the  actual  transition  because  it  applies 
only  while  the  disturbance  is  small  .  The  question  of  how 
'small'  the  disturbance  must  be  in  order  to  apply  the  line- 
arized equations  has  not  been  adequately  answered. 

An  axisymmetric   flow  is  always  stable  to  a  disturbance 
which  has  only  a  tangential  component  of  velocity,  [Ref.  7], 
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Therefore  the  theoretical  analyses  have  generally  been 
limited  to  axi symmet r ical  disturbances  .  Such  analyses  have 
been  carried  out  by  many  investigators,  [Refs.  2,3  and  7]  , 
altough  some  doubt  still  remains  concerning  the  details  of 
the  mathematics. 

Steady  flow  in  a  tube  has  been  found  to  be  stable  to 
small,  axisymmetr ical  disturbances  by  a  number  of  investiga- 
tors,[Refs.  5,7  and  8],  and  may  be  considered  as  an  accepted 
conclus  ion. 

Several  investigators  have  attempted  to  derive  a  cri- 
tical Reynolds  number  ,  for  a  flow  subjected  to  external 
disturbances  ,  after  which  transition  to  turbulence  occurs. 
Leite,  [Ref.  9],  studied  the  behavior  of  small,  axisymmetri- 
cal  disturbances  in  steady  flow  in  a  tube  .  He  found  that 
they  were  always  damped  for  Reynolds  number  up  to  13000  ,  as 
long  as  their  amplitude  was  not  too  large.  Dar 1 ing , [Ref . 10 ] , 
working  with  pulsating  flows  in  tubes  ,  found  that  the  tran- 
sition Reynolds  number  dropped  from  2510.  in  steady  flow  to 
1500  in  pulsating  flow.  Sarpkaya , [Ref .  11],  determined  expe- 
rimentally the  critical  Reynolds  number  for  pulsating  Poi- 
seuille  flow  as  a  function  of  a  frequency  parameter  and 
velocity  ratio  .  Gilbrech  and  Combs,  [Ref.  12],  using  water 
as  working  fluid,  found  a  critical  Reynolds  number  which  was 
higher   than  2220  ,  under  some  conditions,  for  steady  flow  . 
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They  also  pointed  out   that  earlier   investigators  had  found 
no  upper  limit  to  the  upper  critical  Reynolds  numbers. 

These  experimental  results  confirm  the  conclusion  that 
steady  ,  fully  developed  flow  in  tube  is  stable  to  small  di- 
sturbances ,  and  the  critical  Reynolds  number  is  a  function 
of  the  frequency  of  external  disturbances  and  the  physical 
properties  of  the  fluid  .  In  spite  of  this  confirmation  , 
transition  to  turbulent  flow  is  invariably  observed  to  occur 
at  sufficiently  high  Reynolds  number  and  it  is  necessary  to 
consider  this  fact. 

The  most  obvious  explanation  is  that  transition  is 
caused  by  finite  disturbances  which  are  not  damped  ,  even 
though  smaller  disturbances  are  .  This  possibility  is  stron- 
gly supported  by  the  observation  made  by  Gilbrech  and  Combs 
[Ref.  13],  that  the  critical  Reynolds  number  decreases  as 
the  magnitude  of  the  disturbance  increases  . 

The  theoretical  analysis  of  finite  disturbances  has 
been  limited  by  mathematical  difficulties  ,  however  ,  Stuart 
[Ref.  14],  has  shown  that  a  flow  may  be  unstable  to  finite 
disturbances  ,  while  it  is  stable  to  infinitesimal  ones.  In 
addition,  Spielberg  ,  [Ref.  15],  has  pointed  out  that  stabi- 
lity to  two  dimensional  disturbances  does  not  necessarily 
imply  stability  to  three  dimensional  disturbances.  This  lat- 
ter result  creates  a  problem  in  comparing  accurately  the 
theoretical  and   experimental  results  . 
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In  the  present  work,  it  was  assumed  that  an  axisymmetri- 
cal   disturbance   was   propagated  without   confirming   the 
non-exi stance  of  three  dimensional  disturbances  . 

Another  explanation  of  the  origin  of  turbulence  in  pipe 
flow  may  be  that  the  transition  occurs  in  the  entrance 
region  of  the  pipe,  in  which  the  velocity  distribution  is  of 
the  boundary-layer  type  .  A  theoretical  analysis  by  Tatsumi, 
[Ref.  2],  indicates  that  the  flow  in  the  entrance  region  is 
unstable  to  small  disturbances.  He  found  a  point  of  maximum 
instability  near  the  entrance  to  the  tube,  where  the  maximum 
Reynolds  number  for  stability  was  19400  .  However  ,  the  cri- 
tical Reynolds  number  increases  rapidly  in  both  directions, 
so  that  one  cannot  conclude  that  this  local  instability  will 
lead  to  transition.  This  is  because  amplification  required 
to  produce  transition  remains  undetermined. 

Sparrow  and  Lin,  [Ref.  16],  pointed  out,  in  their  effort 
to  determine  the  detailed  nature  of  the  flow  development  in 
the  entrance  region,  that  even  for  a  laminar  flow, the  velo- 
city problem  for  the  entrance  region  does  not  yield  an  exact 
solution,  regardless  of  the  shape  of  the  duct  cross  section. 
Their  explanation  for  the  difficulties  in  the  analysis  was 
that  the  essential  nonlinearity  of  the  inertia  terms  which 
appear  in  the  equation  of  motion,  preclude  linear  analysis. 

From  a  more   sophisticated   point  of  view  ,  the  problem 
becomes   especially   complicated  if  we  attempt  to  reveal  the 
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nature  of  motion  in  the   non-linear  range  of  boundary   layer 
instability  and  the  onset  of  turbulence. 

Klebanoff  and  others  ,  [Ref.  17]  ,  who  attacked  the  pro- 
blem in  the  general  form  ,  pointed  out  that  despite  the  suc- 
cess of  the  linearized  theories  in  revealing  the  nature  of 
the  initial  stages  of  boundary-layer  instability  ,  there  re- 
mains a  deep  void  in  the  undestanding  of  the  subsequent  non- 
linear behavior  and  the  actual  breakdown  of  the  laminar 
boundary-layer . 

Miller  and  Fejer,  [Ref.  18]  ,  who  studied  the  transition 
in  Blasius-type  boundary-layer  produced  by  a  free  stream  ha- 
ving an  oscillatory  component  of  velocity  ,  found  that  the 
transition  Reynolds  number  depends  only  on  the  amplitude  of 
of  the  oscillations  ,  and  that  the  dimens ionless  transition 
length  is  a  function  only  of  the  frequency. 

Obremski  and  Fejer  ,  [Ref.  19]  ,  continued  the  previous 
work  and  pointed  out  that  below  a  critical  value  of  the  pa- 
rameter [  (Re)N<,=  LAU/2rrv  ]  ,  transition  occurs  at  a  relati- 
vely constant  Reynolds  number  ,  which  appeared  independent 
of  the  amplitude  and  frequency  of  the  oscillation  ,  at  least 
over  the  range  investigated. 

Finally  Landahl  ,  [Ref.  25],  used  kinematic  wave  theory 
to  determine  under  what  conditions  breakdown  of  a  steady  or 
unsteady   laminar   flow   into   high   frequency   oscillations 
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should  occur  .  Isolating  the  three  different   ingredients  of 
the  mechanism  of  breakdown  ,  he  pointed  out: 

a.  The   primary   instability  process   its.elf  is  only 
incidental  to  breakdown  and  transition  ,  since  one  can 
easily  produce  a  breakdown  condition   by  a  sufficient- 
ly intense  local   disturbance  of  the   shear  flow  ,  for 
example  through  free  stream  disturbances. 

b.  If  the   breakdown   is  to  be   self  maintained  as  it 
travels  downstream  however  ,  hydrodynamic  instability 
is  required.  In  a  sense  ,  the  breakdown  mechanism  is 
the  one  most  essential   in   transition   of  a  boundary- 
layer  to  turbulence  rather  than  the  classical  hydro- 
dynamic  instability  one,  since  the  former  represents 

a  strongly  irreversible  process. 

c.  An   unstable   small   amplitude   wave  packet  of  the 
Tolmi en-Schl icht ing  type  may  amplify  for  a  while  ,  as 
it  moves  downstream  ,  but  once  it  has  passed  through 
some   streamwise  position  in  the   boundary-layer   it 
will  leave  the  shear  flow  practically  as  undisturbed 
as  before  the  passage  of  the  packet  ,  except   for  a 
small   deformation   of  the  mean  flow  of  second   order 
in  disturbance   amplitude  .  At  breakdown  on  the  cont- 
rary there  is  an   irreversible  redistribution  of  the 
basic  shear  flow  vorticity  due  to  a   nonlinear  recti- 
fication mechanism  . 

In   view  of  the  difficulties   encountered   in  applying 
the  small  disturbance  theory  to  steady  axi symmet r ical  flow  , 
it  is   likely  that   considerable  time  and  effort  will  be  re- 
quired  to   obtain  an  exact   solution  for  unsteady  flow  pro- 
blems ,  especially  if  we  consider  a  three  dimensional  flow  . 

The  use  of  a  '  quasi-steady  *  analysis  would  reduce  the 
problem  to  one  which  could  be  handled  more  easily.  Using 
this  technique  one  would  analyze  an  instantaneous  velocity 
profile  as  though  it  were  a  steady-flow  problem. 
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By  considering  the  stability   of  the  velocity  profiles 
at  various  times,  one  might  then  draw  some  conclusions  about 
the  stability  of  the  unsteady  flow  an  approach  used  by  Shen, 
[Ref.  21],  as  well  as  Greenspan  and  Benney,  [Ref.  22]. 

In  the  present  work,  although  Reynolds  Numbers  as  high 
as  19400  were  not  reached  ,  we  recognize  that  the  flow  in 
the  entrance  region  has  ,  in  all  probability  ,  affected  the 
downstream  flow  and  may  well  have  introduced  turbulence  sub- 
sequently damped  along  the  tube  ,  for  low  Reynolds  numbers  , 
creating  a  laminar  or  laminar-like  flow  at  the  tube  exit. 

B.   NONSTEADY  FLOW  IN  A  CIRCULAR  TUBE 

In  order  to  study  how  the  pressure  disturbances,  intro- 
duced at  the  entrance  of  a  long  tube  used  to  form  the  jet  , 
affect  the  flow  at  the  exit  of  the  tube  ,  which  constitutes 
the  initial  conditions  of  the  flow  of  interest  ,  we  may  exa- 
mine the  nonsteady  flow  in  a  long  tube  of  circular  cross 
section  .  We  adopt,  as  shown  in  Fig.  4  ,  cylindrical  coordi- 
nates whose  x-axis  is  identified  with  the  center  line  of  the 
pipe. 

The  equation  of  continuity  is: 

V(pV)  +  —  =  o 
In  cylindrical  coordinates   Eq.  1   becomes: 
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1  *  1  a  a  "fcp     - 

-—  (rpv  )  +  -  —  (pve)  +  -  (pvx)  +  —  =  0 
p  dp  P  ^9  <>x  0t 


(Eq.   2) 


X ~ 


Figure  4.   Cylindrical  Coordinates  in  a  Round  Pipe. 


Assuming  incompress ibi 1 i ty ,  (o  =  constant),  and  a  non- 
twisting  flow   (v   =  0),   Eq.  2   may  be  written  as: 


1  "a  ,        * 

_    (rpv  )  +     (p   )  =  0 

P  ^P 


d 


or 


c>  v    "^v    V 
x     r     r 

t)  x    c)r    r 


(Eq.   3) 


Neglecting  external  forces  ,  (i.e  gravity)  ,  the  Navier- 
Stokes  equations  in  two  components  can  be  written: 


^x  +    <^x     ^x 
*t     ^x    Vr^r 


1  J)P   p 

p  dx    p 


'bS   ^vx 

Ox2   +  dr2 


1   v. 


+  — 


p    rj 


(Eq.   4) 
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If  we  assume   quasi-steady   flow,  the  equation  of  conti- 
nuity  Eq.  3   becomes: 


it 


=  0 


or:    vx  =  v(r,t) 


(Eq.   6) 


This   indicates   that  the   velocity  in  the  direction  pa- 
rallel to  the  center  line  is  a  constant  at  each  moment. 
Substituting  v   =  0  into   Eq.  5   we  get: 


=0      or:     P  =  P(x,t) 

Ar 


(Eq.   7) 


Inserting   these   results   in  the   equation   of  motion 
(Eq.  4)  : 


p     [^ vx    l      vxl 
x     br2     r   r 


¥*  - . 

1      P 
+   V 

P       x 

v       X 

3t 

h>r2 

where :      V  = 

¥ 

P 

(Eq.   8) 


is  the  kinematic  viscocity  of  the  fluid. 
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C.   SOLUTION  FOR  A  UNIFORM  PULSATING  FLOW 

Assuming  that  the  pressure  pulses  introduced  by  the 
drive  unit,  (a  loudspeaker  driver),  are  propagated  instantly 
along  the  flow  direction  ,  (  justified  by  the  fact  that  the 
mean  flow  velocity  is  very  small  compared  with  the  speed  of 
sound  in  the  fluid  )  ,  with  the  aid  of  Eqs.  6,  7  and  8  we 
can  show  that  the  pressure  gradient  ,  (  c>  P/  3  x  ),  becomes  a 
function  only  of  time  ,  t. 

Using  Fourier   series  we  can  express  the  pressure  gra- 
dient as  follows  : 

1  }P  ~  ^ 

=   or     +   {&    ocncos(nt)  +    ^j    bnsin(ht)  (Eq.      9) 

p    clX 

or,  for  the  sake  of  simplicity  of  calculation  ,  in  a  complex 
form: 

up  -5-     •  t 

where  :  oc  ,  b   are  constants  representing  the  amplitudes  of 
elemental  vibrations. 

The  corresponding   solution  for  the  axial   velocity  may 
be  assumed  to  be  of  the  form  : 


vx    =    V„    +    ^     v„    nosfntf   +    ^4     v_     «  i  nfn  t\  (Eq.     11) 


vx0+    n=\    vnicos(nt)+    ^    vn2sin(nt) 
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or,  in  complex  form: 

so 

vx  -  vXo+  n?!  vneint   .   vn  -.  ^  -  u^         (Eq.  12) 

where  coefficients   v   and   v   are  assumed  to  be  functions 

n4       nz 

of  the  radial  coordinate  ,  r  ,  only. 

Substituting  Eqs .  10  and  12  into  Eq.  8  and  comparing 
the  terms  of  the  same  order  ,  we  obtain  the  following  Dif- 
ferential Equations  for  the  coefficients  as  independant 
var  iables. 


d2v    i  dvx   vXo 
dr2     r  dr     V 


and : 


After  integration  we  get: 

2 
4 


Xo 


where:    A   ,B   =  constant 

o  '  o 

and  : 


v 

l  n 


(Eq.  13) 


d2vn    1    dvn    "in     ^n 

—2  + +  vn  + =  0  (Eq.  14) 

dr2    r  dr     V      v 


ex.     p 
= -  —  +  Aoin(r)  +  B0  (Eq.  15) 


n  ■  T-  +  DnJ0(o-ri3/2)  +  EnK0(5-ri1/2)  (Eq.  16) 
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where:   5  =']/  (n/v) 

and  Jo  and  Ko  are  the   Bessel   functions  of   the   first   and 

second  kind  of  zero  order  respectively. 

The  solution  for  the   x-wise   component  of  velocity  is 
given  by  the  real  part  of  the  following  expression: 


ao    r 

vx    = +  AQlnr    +   B0   + 

x  v      4  °  ° 


oo 


n  =  l 


fi  +  Dnj0(§:ri3/2)  +  EnKo(S:rii/2) 


eint  (Eq.    17) 


On  the  axis  of  the  pipe  ,  (  r  =  0  )  ,  ln(r)  =  -oo   and 

1/2 

Ko(8tri         )    =  oo  ,    since      vx   is    finite    ,   we   conclude    that: 

A      =    E      =    0  (Eq.    18) 

on 

The  Boundary  Condition  at  the  wall  is  given  by: 

v   =0   at   r  =  R  (Eq.  19) 

where:   R  is  the  radius  of  the  pipe. 

Inserting   this   boundary   condition   ,  Eq.  19  ,  into 

Eq.  17  we  get: 


B   = 
o 


(XoRi 
4v 


n=l 


Lin 


(Eq.  20) 


Inserting  ,  Eqs .  18  and  20  into  Eq.  17,  we  finally  get 
the  velocity  component  v   : 
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or 

i 

4v 


oo 
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1    - 


J0(Sri3/2) 

J    (5Ri3/2)l 
o 


int 


(Eq.    21) 


The  total  mean  mass  flow  ,  m  ,  which  is  transmitted  in 
the  x-direction  is  given  by  the  integral  over  time  and  ra- 
dius of  the  axial  velocity  component: 


!      2  tt        R 
m  =  /dt/2  Trvxrd 


TrR4a 


r  = 


8v 


(Eq.  22) 


Inserting  the  time  mean  pressure  gradient: 


bx 


=    Pao 


(Eq.    23) 


we    get : 


m   = 


MR 


8u 


4, 


<>P 

Dx 


(Eq.    24) 


It  is  seen  that  the  total  mean  mass  flow  ,  in  pulsating 
motion  %    is  identical  to  that  of  a  steady  Poisseuille  flow 
with  the  same  pressure  gradient,  as  the  mean  pressure  gra- 
dient in  the  pulsating  flow. 

A  mean  velocity  ,  U  ,  may  be  defined  by: 


U  = 


m 


TTR' 


(Eq.  25) 
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Combining  Eqs.  24  and  25  we  get 


R2  f  dP 
U  = -  — 

8p  I   Ox 


(Eq.  26) 


We  may  express  the   x-wise  velocity   as  the  sum  of  a 
perturbation  velocity,  v^  f  an(j  tne   mean  stream  velocity  of 

the  fluid,vxs  : 


vx  =  vxs  +  vx 


(Eq.  27) 


Introducing  the  Hagen-Poi seui 1 le  relation  for   v    we 

A  S 


get: 


xs 
U 


=  2 


r2! 
1 


(Eq.  28) 


Dividing   Eq.  27   by   the  total  mean  velocity  ,  U  ,  we 
get  the  non-dimensional  expression: 


v    v     v 
x  _   xs     x 

U     U    U 


(Eq.  29) 


Combining      Eqs.    21    and    29   we    get    for    the    perturbation 


velocity    ,    v      : 
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*  OO 


anr    8B 
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n=l 


OO 
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8B 


(&R) 


8(1-A) 

-cosfrit)  +  9    sinfrit) 

z  ( Sr ) L 

8(1-A) 
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2  ( SR ) L  J 


(Eq.    30) 
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where: 


A  = 


B  = 


ber(§R)ber(o"r  )  +  be i (o"R)be i (£r ) 

ber2(o*R)  +  bei2(o~R) 

bei  (<5R)ber  (Sr)  -  ber (o*R)be i (Sr ) 

ber2(5R)  +  bei2(5k) 


(Eq.  31) 


where:   ber   and   bei  are  the  Kelvin  functions  . 

With  the  aid  of   Eqs.  10  and  26  we  can  obtain  the  cor- 
responding pressure  gradient  relation: 


4R  ^P    64 


OO 


©o 


pU  ^x   Re 


1  +  -^  —cos(nt)  +  ±T,   — sin(nt) 
n=1  a         n_1  o< 


(Eq.  32) 


where: 


Re  =  2RU  /  V 

Introducing  the  friction  factor  for  laminar  flow  : 

f  =  64  /  Re 

into  Eq.  32,  we  get  : 


^P        pU2 


4R 


1    +   n^i—  cosfrit)  +   ^i  —  sinfrit) 


(Eq.  33) 


where:  an/o;0  ancj  bn/ofo  are  dimens  ionless  ratios  of 
amplitude  of  periodic  variation  of  pressure  gradient  norma- 
lized with  the  mean  amplitude  . 
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D.   VELOCITY  DISTRIBUTION  AND  ASYMPTOTIC  EXPRESSIONS 

The  velocity  profiles  for  laminar  ,  pulsating  flow  cons- 
titute the  initial  conditions  for  any  theoretical  analysis 
of  stability.  Although  no  exact  theoretical  analysis  of  sta- 
bility is  presnted  in  this  work,  solution  for  the  velocity 
profiles  will  serve  to  identify  the  dimens ionless  parameters 
which  are  significant  and  will  permit  some  qualitative  de- 
duction concerning  the  stability  of  pulsating  flow. 

From  Eq.  16  we  found  that  the  periodic  part  of  the  velo- 


city profile  is  characterized  by  the  parameter   Sr  =|/n/v  R. 

The  distribution  of  ocn/  <x0    ,  the  dimens  ionless  velocity 
amplitude,  obtained  by  Gilbrech  and  Combs  ,  [Ref.  16]  ,  is 
plotted  in  Figure  5. 

Uchida,  [Ref. 23],  obtained  solutions  for  limiting  values 
of  the  parameter  oR  by  means  of  asymptotic  expansions  of  the 
Bessel  function: 


1.   Slow  Oscillation  (  SR  <<  1  ) 

When  a  highly  viscous  fluid  oscillates  slowly  in  a 
narrow  pipe  ,  oR  becomes  small  and  we  can  approximate: 


ber(oR)  —  1    and   bei(SR)  —  0 

Using  this  approximation,  we  get  from  Eqs.  28  and  29 
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Figure  5.   Oscillation  Amplitude  of  the  Local  Velocit 
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or : 


=  2 
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R2  J  ot  I 


1  ^P 

p  "dx. 


or : 


1  o   f   X    ^Pl 

s_(R»-r»>  -  — 
4V         I   P  ^xj 


(Eq.  34) 


1  ^P 

4V         l   p 

From  this  we  conclude  that  the  velocity  distribution 
is  a  parabola,  as  in  the  case  of  steady  Hagen  -  Poisseuille 
flow,  while  the  magnitude  of  velocity  varies  periodically  in 
phase  with  that  of  the  pressure  gradient. 


2.   Rapid  Oscillation   (  5R  -»-  °°  ) 

When  a  fluid  of  low  viscocity  oscillates  rapidly  in 
a  large  pipe  ,  then  parameter   oR  becomes  large. 

Assuming  5R  >  10  ,  asymptotic  expansion  of  Bessel 
functions  may  be  introduced: 
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ber(SR)  = —  *fcos 
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>  (Eq.  35) 
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To  examine  the  motion  of  the  fluid  near  the  center 
of  the  pipe  ,  we  put: 

Sr  — .  oo   and   Sr  — —  0 


As  before  we  get: 
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or : 
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(Eq.    36) 

From  this  result  we  may  conclude  that  in  the  case  of 
rapid  oscillation  ,  fluid  flows  at  the  center  of  pipe  with 
the  phase  lag  of   90    behind  the  wave  of  pressure  gradient 
and  its  amplitude  diminishes  with  increasing  frequency. 

To  examine  motion  of  the  fluid  near  the  wall  of  the 
pipe  ,  we  put: 

5R  -*-o©    and    6>R  -*eo 

Using  the  same  procedure  we  get: 
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^2 


Because  of  the  obvious  complexity  of  the  expression, 
it  is  difficult  to  draw  any  simple  conclusions  other  than  to 
note  that  a  local  maximum  in  velocity  developes  in  the  neigh 
borhood  of  the  wall  as  it  is  shown  in  Fig.  5  .  We  can  also 
get  a  better  confirmation  from  Figs.  6  ,  7  and  8  ,  obtained 
from  [Ref.  23],  for  values  of  5R  =  1 ,  3  and  10  respectively. 


39 


0 

"a 


n/ 


0°     - 

_ 

>^~ 

^N 

330°"" 

JX.-1 


Figure  6.   Velocity  Profile  for  $R   =  1 


40 


Figure  7.   Velocity  Profile  for   5R  =  3 
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Figure  8.   Velocity  profile  for   Sr  =  10 
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E.   SECTIONAL  MEAN  VELOCITY 

Since  the  instantaneous  mass  flow  and  accordingly  sec- 
tional mean  velocity  ,  ape  periodic  functions  ,  we  seek  to 
link  the  mean  velocity  with  the  applied  pressure  gradient. 

Sectional  mean  velocity  ,  denoted  by  vXJrt,  is  given  by: 

1      } 

727TVxrdr  (Eq.    38) 


_TtR2^27TVxrdr 
o 


Introducing  Eq.  34  and  integrating  we  get: 
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(Eq.    39) 


d[ber(o"R)l               ,                d[ber(o~R)] 
ber'(SR)    =-i ,    bei(oR)    =   (Eq.    40) 
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(Eq.41) 


An  example  of  simple  periodic  pulsation  is  given  by: 
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(Eq.    42) 


where : 
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(  average  term  ) 


=  a  cosnt     (  perturbation  term  ) 


Inserting  the  above  relation  into  Eq.  39  we  get: 
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(Eq.    43) 
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is  the  coefficient  of  amplitude  ,  and 
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represents  the  coefficient   of   phase   lag   in  the  wave  of 
pressure  gradient. 

In   Fig.  9  ,  obtained  from  [Ref.  23]  ,  are  plotted  the 
coefficients  of  amplitude  and  phase  lag  of  mean  velocity  as 
a  function  of  the  parameter  oR. 
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Figure   9.   Coefficients  of  Amplitude  and  Phase 
Lag  of  Mean  Velocity. 

F.   SURFACE  FRICTION  AND  BALANCE  OF  FORCE 

The  instantaneous   frictional   force  acting  on  the  pipe 
wall  is  given  by: 


T  =  -1 


dv 


dr 


r=R 


(Eq.  44) 


Dividing   Eq.  44   by  (  pU  /2  )  ,  and  combining  with  the 
previously  obtained   Eqs .  32  ,  39  and  40  we  can  get: 
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(Eq.  45) 


Considering  the  simple  periodic  pulsation  described  by 
Eq.  42  we  get : 
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(Eq.  46) 


represents  the   coefficient  of  amplitude  of  shearing  stress 

at  the  wall  ,  and: 


$T  =  tan-1 


D 

L  C 


represents  the   coefficient  of  phase  lag  of  shearing  stress 

with  respect  to  the  wave  of  pressure  gradient. 
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In  the  Fig.  10  below  ,  obtained  from  [Ref.  23],  it  is 
shown  how  the  coefficient  of  amplitude  and  phase  lag  of 
shear  stress  varies  with  %R 


Figure  10.   Coefficients  of  Amplitude  and  Phase  Lag 
of  the  Shearing  Stress 

If  we  compare  Figs.  9  and  10  we  observe  no  substantial 
differences  in  amplitude  or  phase  lags.  Moreover,  the  phase 
of  the  sectional  mean  velocity  is  greatly  delayed  with  res- 
pect to  the  pulsating  wave  of  pressure  gradient,  while  that 
of  shearing  stress  is  less  delayed  from  it. 

Using   Eq.  8   and   integrating   over  the  section  of  the 
pipe,  from  r=0  to  r=R  ,  we  can  obtain  a  relationship  between 
the  force  and  acceleration  of  amass  of  fluid  enclosed  in  the 
circular  cylinder  of  unit  length: 

J>P  dv 

-  TTR2 =     TTR2p — —  +    2-ttRt 

c)x  dt 


47 


1    -  dP/dx 

or :     — 


dvxm/dt 


4r 


2      pD2/4R  D2/2R 


(Eq.  47) 


PDi 


For  the  present  example  of  a  simple  periodic  pulsation 
we  find  the  mean  frictional  force  ,  involved  in  the  shear 
stress  ,  is  just  balanced  by  the  term  of  mean  pressure  gra- 
dient and  there  results  no  steady  acceleration.  The  various 
terms  appearing  in  Eq.  47  are  shown  in  Fig.  11  below  ,  ob- 
tained from  [Ref.  23]. 
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Figure  11.   Periodic  Force  Component  (A  -  Pressure  Gradient, 
B  -  Acceleration  ,  C  -  Shearing  Force) 
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IV.   RESULTS  AND  CONCLUSIONS 

A.   PRELIMINARY  RESULTS 

After  a  preliminary  survey  of  the  phenomenon  of  jet 
breakdown  at  a  single  laminar  and  single  turbulent  Reynolds 
number,  as  mentioned  above,  it  was  found  that  the  field  was 
so  rich  in  information  that  a  third  ,  transitional  Reynolds 
number  of  2200  was  added  to  the  survey  .  Observations  repor- 
ted are  of  the  most  dinstinctive  phenomena  at  each  of  the 
three  Reynolds  numbers  and  were  observed  to  occur  most  fre- 
quentlyat  frequencies  which  were  multipliles  of  33  and  50  hz 
Results  include  two  photographs  at  each  data  point  of  two 
separate  phases  along  with  an  oscillogram  showing  both  the 
excitation  waveform  and  the  pressure  perturbation  it  produ- 
ced, shown  in  Figs.  13  through  36.  Table  I  is  a  summary  of 
the  observed  breakdown  length  as  a  function  of  the  frequency 
and  Table  II  is  a  summary  of  oscilloscope  and  flow  parame- 
ters, corresponding  to  each  figure. 

As  Hoyt  and  Taylor,  [Ref.  24],  pointed  out,  jets  are  un- 
stable to  two  kinds  of  disturbances,  axial  and  helical.  The- 
ory suggests  that  the  axial  instability  (  corresponding  to 
the  initial  waves  on  the  jet  surface  )  die  out  after  a  few 
nozzle  diameters  downstream  and  the  helical  instabilities, 
amplified   by  aerodynamic  resistance  ,  lead  to  the  final  jet 
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breakup.  The  initial  slight  helical  motion  is  amplified  un- 
til the  entire  jet  may  be  described  as  having  a  corksrew- 
like  motion  .  This  is  well  confirmed  by  the  photos  obtained 
in  the  present  work  of  the  jet  formation  before  and  after 
the  breakdown  and  is  in  good  agreement  with  the  observations 
obtained  by  Crow  and  Champagne,  [Ref.  4], 

TABLE  I 


SUMMARY  OF  BREAKDOWN  DISTANCE 


f  (Hz) 

#d 

#d 

#d 

(Re=1760) 

(Re=2200) 

(Re=2640) 

0 

118.0 

125.0 

30.0 

33 

112.6 

120.6 

31.5 

50 

110.0 

112.0 

34.0 

66 

106.3 

103.7 

32.2 

100 

102.0 

98.5 

29.5 

133 

98.6 

99.2 

33.7 

150 

95.5 

100.0 

35.0 

166 

92.7 

101.7 

37.8 

200 

88.0 

102.5 

41.0 

233 

85.2 

101.2 

45.2 

250 

83.5 

99.0 

49.0 

266 

81.5 

101.5 

47.6 

300 

91.2 

88.0 

45.0 

333 

96.1 

97.5 

44.1 

350 

101.8 

110.3 

43.3 

366 

110.0 

121.0 

42.0 

400 

112.0 

125.0 

42.5 

Re  =  Reynolds  Number 

f   =  Introduced  Frequency 

#d  =  Number  of  Diameters 
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TABLE  II 


SUMMARY  OF  EXPERIMENTAL  CONDITIONS 


Fig.# 

Re 

f 

V 

\r 

P 

Vg 

t /cm 

(Hz) 

(mV/cm) 

(mV/cm) 

(msec/cm) 

13 

1760 

50 

5 

10 

5 

14 

1760 

100 

5 

10 

5 

15 

1760 

150 

5 

10 

2 

16 

1760 

200 

5 

10 

2 

17 

1760 

250 

5 

10 

2 

18 

1760 

270 

5 

10 

2 

19 

1760 

300 

5 

10 

2 

20 

1760 

380 

5 

10 

2 

21 

2200 

50 

5 

10 

5 

22 

2200 

100 

5 

10 

5 

23 

2200 

150 

5 

10 

2 

24 

2200 

200 

5 

10 

2 

25 

2200 

233 

5 

10 

2 

2  6 

2200 

250 

5 

10 

2 

27 

2200 

266 

5 

10 

2 

28 

2200 

300 

5 

10 

2 

29 

2200 

380 

10 

10 

1 

30 

2200 

787 

20 

20 

0.5 

31 

2640 

50 

5 

10 

5 

32 

2640 

100 

5 

10 

2 

33 

2640 

200 

5 

10 

2 

34 

2640 

250 

10 

10 

2 

35 

2640 

300 

10 

10 

1 

36 

2640 

600 

10 

10 

0.5 

Re  :  Reynolds  number 

f  :  Introduced  Frequency 

p  :  Signal  Amplitude  from  Pressure  Transducer 

V 

g  :  Signal  Amplitude  from  Frequency  Generator 
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Figure  12.   Breakdown  Length  Versus  Frequency 
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1.   Laminar  Flow  at  Reynolds  Number  1760 

At  this   Reynolds  number   the  unexcited  flow  started 
to   breakdown  at  approximately   118   diameters  from  the  tube 
exit.  Increasing  the  frequency  to  266  Hz  produced  a  gradual- 
ly movement  of  the  breakdown  point  upstream  to  approximately 
81,5  diameters  from  the  exit.  Above  this   frequency  a  rever- 
sed trend   appeared   and  the  breakdown  point  started  to  move 
downstream  as  the  flow  became  laminar  similar  to  the  unexci- 
ted flow  ,  the  only  difference  being  in  the  droplets  formed 
after  the  breakdown.  This  phenomenon  may  be  characterized  as 
a   re laminar izat ion   of  the  flow  and  was  almost  complete  at 
380  Hz.  Further   increasing  of  the  frequency  up  to   1000  Hz 
produced  no  significant  effect  in  the  jet  flow. 

Results  for  Reynolds  number  1760  are  shown  in  Figs. 
13  through  20  .  As  can  be  observed  from  the  photos  ,  break- 
down point  is  not  fixed  and  oscillates  somewhat  upstream  and 
downstream  ,  creating  a  breakdown  region  .  The  amplitude  of 
this  oscillation  was  reduced  as  the  perturbation  frequency 
was  increased  ,  starting  from  about  14  diameters  and  approa- 
approching  4  diameters  as  the  frequency  changed  from  200  to 
266  Hz  .  At  266  Hz  an  especially  interesting  phenomenon  was 
observed  concerning  the  breakdown  point  and  the  droplet  for- 
mation .  The  breakdown  of  the  jet  occured  at  the  minimum  ob- 
served distance  from  the  tube  exit  at   81.5   diameters. 
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Figure    13.      Photos    Illustrating   Jet    Breakdown    and     iscilloscope 
Trace   of    Pressure   at    Re=1760   and    f=50  Hz 
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Figure    14. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=1760  and  f=100  Hz 
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Figure  15.   Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=1760  and  f=150  Hz 
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Figure    16. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=1760  and  f=200  Hz 
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Figure  17. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=1760  and  f=250  Hz 
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: 


Figure  18. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=1760  and  f=266  Hz 
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Figure  19.   Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=1760  and  f=300  Hz 
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Figure  20. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=1760  and  f=380  Hz 
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The  oscillation  of  the  breakdown  point  was  minimized  to  a- 
bout  4  diameters  and  the  droplets  formed  ,  as  was  observed 
with  the  aid  of  the  stroboscope  ,  were  exactly  the  same  and 
appeared  stationary  in  space.  This  is  shown  in  the  jet  pho- 
tos of  Fig.  18. 

A  plot  of  the   breakdown   distance  as  a  function  of 
frequency  is  provided  in  Fig.  13  where  we  see  that,  in  the 
frequency  range  from  zero  to  266  Hz,  an  almost  linear  varia- 
tion of  the  breakdown  point  with  frequency  occurs.  The  mini- 
mum at  266  Hz  may  be  characterized  as  a  critical  frequency 
for  this  specific  Reynolds  number  .  Above  this  frequency  one 
also   observes   an  almost   linear  variation  of  the   average 
breakdown  point  to  380  Hz.  This  region  may  be  characterized 
as  relaminar izat ion  region.  The  slope  of  the  breakdown  curve 
is  greater  than  at  lower  frequencies  although  increasing  the 
frequency   above   380  Hz   results   a  zero  slope  curve  ,  as 
discussed  above. 

2 .   Transitional  Flow  at  Reynolds  Number  2200 

At  this  Reynolds  number  the  unexcited  flow  started 
to  breakdown  at  approximately  125  diameters  from  the  tube 
exit  .  Increasing  the  frequency  to  100  Hz  ,  the  breakdown 
point  moved  upstream  up  to  about  98  diameters  from  the 
tube  exit  .  The  oscillation  of  the  breakdown  point  was 
reduced  at  this  case,  from  12  to  6  diameters. 
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Increasing   the  frequency   to   266  Hz   resulted  in 
only   slight   variation  of  the   breakdown   point  .  The  only 
unusual   observation   in  this   range   of  frequency  was  the 
drastic   reduction   of  the   amplitude   of   breakdown   point 
oscillation  as  the  frequency  was  increased  . 

An  interesting  pressure  effect  was  observed  at  266 
Hz,  as  can  be  seen  at  Fig.  27  ,  where  a  maximum  amplitude 
of  pressure  signal  was  recorded. 

An   interesting   change  was  observed   when  the  fre- 
quency reached   300  Hz  .  The   breakdown  point  was  displaced 
suddenly  about  25  diameters  upstream  and   nearly  stabilized 
at  88  diameters  .  The  droplets   viewed   with  the  aid  of  the 
stroboscope   appeared   immobile   and  any  oscillation  of  the 
breakdown  point  was  difficult  to  observe  .  Based  on  a  series 
of  high   speed   photos  it  was   concluded   that  the  breakdown 
point  was  oscillating  with  an  amplitude  of  about  2  diameters 
Also  observed  at  this  frequency  that  the   waveforms  obtained 
from  the  signal  generator   and  the  pressure  transducer   were 
almost  in  phase  and  moreover,  the  latter   waveform   appeared 
nearly  sinusoidal  as  shown  in  Fig.  28. 

Further  frequency  increases  resulted  in  a  flow  rela- 
minarization  as  was  observed   and  discussed  in  the  previous 
case.  A  completely  laminar  flow  was  observed  at  a  frequency 
of  380  Hz  and  no  additional  significant  effects  were  obser- 
ved up  to  2000  Hz  . 
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Figure  21.   Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2200  and  f=50  Hz 
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Figure  22. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2200  and  f=100  Hz 
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Figure  23.   Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2200  and  f=150  Hz 


66 


Figure  24. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2200  and  f=200  Hz 
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Figure    25. 


Photos    Illustrating   Jet    Breakdown    and  Oscilloscope 
Trace   of    Pressure    at   Re=2200   and    f=233  Hz 
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Figure  26. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2200  and  f=250  Hz 
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Figure    27.      Photos    Illustrating    Jet    Breakdown    and  Oscilloscope 
Trace   of    Pressure    at   Re=2200   and    f=266  Hz 
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Figure  28. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2200  and  f=300  Hz 
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Figure  29. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2200  and  f=380  Hz 
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Figure  30.   Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2200  and  f=787  Hz 
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Refering  to   Fig.  12  ,  we  may  draw  some  conclusions 
concerning  the  average  breakdown  point.  Up  to  100  Hz  a  gra- 
dually movement  of  the  breakdown  point  upstream  is  observed 
which  is  nevertheless  more  rapid  than  the  laminar   case  as 
can   be  seen   from  comparison   of  the   slope  of  the  curves 
within  this  region  .  A  nearly  stable  region  then  appears  in 
the   frequency   range  between   100  Hz  and  266  Hz   where  the 
breakdown  point  is  nearly  fixed  at  100   diameters.  This  re- 
gion may  be  characterized   as  a   stabilized  region   and  was 
not  observed  at  any  other  Reynolds  number.  A   singular  point 
was  also  observed  at   300  Hz   as  discussed  previously.  This 
frequency  may  be  characterized  as  a   critical  frequency   for 
the  specific   Reynolds  number  .  Finally  ,  we  can   see   from 
Fig,  12  that  the  relaminar i zat ion   rate  was  also  greater  for 
transitional  flow,  producing  a  steeper  slope  of  the  curve  in 
the  frequency  range  from  300  Hz  to  400  Hz. 

3.   Turbulent  Flow  at  Reynolds  Number  2640 


At  this  Reynolds  number  it  was  observed  that  the 
unexcited  flow  began  to  breakdown  at  an  average  distance  of 
30  diameters  and  no  significant  effect  on  the  jet  was  obser- 
rved  up  to  a  frequency  of  100  Hz.  Increasing  the  frequency  a 
gradual  movement  of  the  average  breakdown  point  downstream 
observed,  which  reached  a  maximum  of  50  diameters  at  250  Hz. 
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Further   increase  of  the  frequency   produced  a  very 
slow  movement  of  the  average   breakdown  point  upstream  to 
about  42   diameters  at  a  frequency  of  366  Hz.  The  breakdown 
point  then   remained   almost  unchanged  up  to   2000  Hz  .  The 
only  further   change   observed  in  the  high  frequency  range  , 
from   366  Hz  to   2000  Hz,  was  the  appearance  of  a  corksrew- 
like  motion  of  the   jet  which  became  more   pronounced  with 
increasing  frequency  ,  making  more  readily  observable  the 
helical   mode  ,  which   dominates  the   breakdown  mechanism  . 
These   results  are  shown  in  Figs.  31  through  36. 

B.  CONCLUSIONS. 

From  the  results  ,  discussed  above  ,  we  may  draw  the 
following  conclusions  concerning  the  effect  of  small  pertur- 
bations on  jet  breakdown: 

a.  In  the  case  of  laminar  flow,  (Re=1640)  ,  the  breakdown 
point  moves  upstream  as  the  frequency  is  increased  up 
to  266  Hz.  Downstream  of  the  breakdown  point  the  jet 
flow  becomes  turbulent  and  breaks  into  discrete  dro- 
plets. Above  266  Hz  this  effect  is  observed  to  reverse 
itself  and  the  jet  returns  to  its  laminar  nature. 

b.  In  the  case  of  transitional  flow  ,  (Re=2200)  ,  the 
breakdown  point  moves  initially  upstream  as  frequency 
is  increased  to  100  Hz  .  In  the  range  of  100  Hz  to 
266  Hz  an  almost  stable  behavior  of  the  breakdown 
point  is  observed  .  Singular  behavior  with  very  inte- 
resting results  including  a  sudden  shift  of  breakdown 
point  is  observed  at  300  Hz  .  Finally  a  more  effective 
re  laminar i zat ion  than  observed  at  lower  Reynolds  num- 
ber occurs  above  300  Hz  leading  to  a  completely  stabi- 
lized laminar  flow  at  frequencies  above  380  Hz  . 

c.  In  the  case  of  turbulent  flow  ,  (Re=2640)  ,  no  notice- 
able  effect  is   observed  up  to  a  frequency  of  100  Hz. 
Above   this   frequency  the  breakdown  point   begins  to 
move  gradually   downstream  up  to  250  Hz  ,  returns  back 
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upstream  to  a  frequency  of  333  Hz  ,  and  becomes  almost 
completely  unaffected  by  the  pressure  perturbations 
at  higher  frequencies. 

d.  Changing  the  waveform  of  the  disturbance  produced 
only  a  slight  deviation  in  the  average  breakdown 
point,  however  a  markedly  different  patern  of  droplet 
formation  was  observed. 

e.  Although  it  was  not  possible  to  record  the  phenomena 
photographically  for  inclusion  here  ,  it  was  possible 
using  the  phase  angle  control  incorporated  in  the 
stroboscope  to  observe  the  entire  development  of  the 
droplet  formation  visually.  This  process  was  found  to 
be  a  complex  sequence  of  breakup  and  re-agglomeration. 
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Figure  31. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2640  and  f=50  Hz 
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Figure  32. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2640  and  f=100  Hz 
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Figure  33. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2640  and  f=200  Hz 
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Figure  34. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2640  and  f=250  Hz 
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Figure  35. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2640  and  f=300  Hz 
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Figure  36.   Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2640  and  f=600  Hz 
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APPENDIX.    CALCULATIONS  AND  DATA  REDUCTION 

A.   CALCULATION  OF  CRITICAL  REYNOLDS  NUMBER 

A  visual  observation  of  the  flow  was  used  in  order  to 
determine  the  critical  Reynolds  number  .  Adjusting  the  flow 
by  the  means  of  the  needle  valve  ,  installed  at  the  entrance 
of  the  plenum  chamber  ,  a  laminar  jet  was  observed  up  to  the 
flow  meter  indication  m.  =  25  .  Above  this  value  a  turbulent 
flow  resulted. 

Using  the  equation  of  continuity  : 

m  =  pAU   or   U  =  m  /  pA  (Eq.  A-l) 

where  : 

m  =  the  fluid  mass  rate. 

U  =  the  mean  flow  velocity. 

p  =  62.34  lbm  /ft3  ,  the  water  density  (at  60°  F). 

A  =  the  cross  section  area  of  the  pipe. 
Also: 

A   =    Ttd2    /    4    =    (tt/    4)[(3    /    16)(1    /    12)]2    (ft2) 
or:    A  =    1.9175    .10~4    (ft2) 
and  : 

m  =  (.6(25/100)](gal/min)  or  m  =  .15  gal/min 
or:  m  =  ( . 15 ) ( 500 . 8  )  ( 1/3600  )  (lbm/sec) 

or:  m  =  .020867  (lbm/sec) 

Substituting  the  above  values  into   Eq.  A-l   we  get: 
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.020867 

U  = (ft/sec) 

(62.34X1.9175  10"4 


or:  U  =  1.74562  (ft/sec) 

The  Reynolds  number  is  given  by: 

Re  =  (pUd)  /  (ugc)  (Eq.  A-2) 

where : 

u  =  2.4  .10"5  (lbf-sec/ft2) 


and : 


g   =  32.2  (lbm-ft )/(lbf-sec2) 
°c 

Substituting  the  above  values  into  (Eq.  A-2)  we  get: 
(62.34)(1.74562)(l/64) 


(Re) 


CF       (2.4  10"5) (32.2) 


or:  (Re)    =  2200 
cr 

The  value  obtained  may  be  considered  reasonable  taking 
into  account  that  is  the  most  common   value  found  in  the  Hy- 
drodynamics  Bibliography.   The  numerical  values   used  in 
the  above  calculations  were  obtained  from  Ref.  25  . 

B.   CALCULATION  OF  INDICATED  PRESSURE  IN  LAMINAR  FLOW 

Refering  to   Fig.  B-l,  where  the  exact  dimensions  used 
are  indicated,   we  apply  the  Energy  Equation: 

Pl      Ui2    B       P2      U22    g 
—  +  b  -±-  +  -  z   =  —  +  b2-^-    +   -  z2  +  hl  (Eq.  B-l) 

P       2go    g0  p       2g     gc 
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Figure  B-l.   Sketch  of  the  Apparatus. 


where  : 

P..  :  the  pressure  at  point  1. 

P0  :  the  pressure  at  point  2  (  P0  =  P    =  14.7  psia) 

z  z     8.  i  m 

U1  :  the  mean  velocity  of  the  fluid  at  point  1. 
U,  :  the  mean  velocity  of  the  fluid  at  point  2. 
b1  =  b„  =  2  :  the  energy  correction  factor  for  laminar 
flow  at  points  1  and  2  respectively. 

z..,z?  :  the  height  difference  from  the  reference  plane 
h,  :  the  flow  losses  given  by  the  relation: 

,2 


hl  "  f 


1   U' 


d  2g, 


(Eq.  B-2) 
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where: 

f  =  64/Re  :  friction  factor  for  laminar  flow. 

1  :  the  total  pipe  length. 
Inserting  the  numerical  values  and  the  appropriate  di- 
mensions from  Fig.  B-l  into  Eq.  B-2  we  get: 

hx  =  (64/2200M94. 416/. 1875)  [( 1 .  74562  )  2/2  (  32  .  2  )] 
or:    hj  =  .69335  (lbf-ft/lbm) 

Using  the  Continuity  Equation   refering  to  Detail  A  of 
Fig.  B-l  and  neglecting  the  minor  turning  losses  we  get: 

m  =  pAUx  =  pAU2  (Eq.  B-3) 

or:    U2/Ux  =  Ax/A2  =  d^/d^  =  (  .25  )  2/ (  .  1875  )  2  =  1.7777 
Therefore : 

Ux  =  U2/1.7777  =  1.74562/1.7777  (ft/sec) 
or:    Ux  =  .982  (ft/sec) 

From  the  Energy  Equation,  (Eq.  B-l)  we  may  write: 

-^  (Px  -  P2)  =-  (V'V*  +  ^VZ1}  +  gchl 
P  2 

Solving    for     AP   =   pi"P2      we   Set : 

iP  =  V  (o^  v  8urt>  *«*Hi      (EqB.4) 

Substituting  numerical  values  into  Eq.  B-4  we  get: 

AP    a!i±l/a     [(1. 74562)2. (.982)2]  +  32.2||  +  32.2(. 69335)} 
32.201 2 

or:        AP   =    125.18    (lbf/ft2) 
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or:        AP   =    .8693      (psig) 

or:       AP   =    1.775      (in.Hg)      at    60°   F 

The  latter  is  a  reasonable   result   compared  with  the 

mercury  manometer  which  indicated  pressure   P.  .  =  1.8  in.Hg 

J  v  ind  & 

suggesting  an  error  of  1.39  %    . 

The  curvature  of  the  pipe  does  not  create  a  large  effect 
in  friction  factor  for  this  Reynolds  number  which  is  in 
good  agreement  with  the  experimental  work  of  White  [Ref.26], 

C.   DYNAMIC  RESPONSE  OF  INDICATED  PRESSURE 

Since  the  pressure  perturbations  were  introduced  in 
the  upstream  plenum  chamber  and  measured  with  a  pressure 
transducer  through  a  pressure  tap,  we  may  reasonably  inquire 
into  the  accuracy  of  the  pressure  traces  measured.  In  order 
to  answer  this  question,  we  consider  the  dynamic  response  of 
the  system  under  the  assumption  that  the  fluid  pressure  is 
applied  uniformly  at  the  transducer  front  surface  . 

The  transient   response  of  a  pressure   measuring  instru- 
ment is  dependent  on  two  factors: 

a.  The  response  of  the  transducer  element  that  senses 
ses  the  applied  pressure. 

b.  The   response   of  the  pressure   transmitting  fluid 
contained  in  the  connecting  tubing. 
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This  latter  factor  is  frequently  the  one  which  deter- 
mines the  overall  frequency  response  of  a  pressure  measure- 
ment system  ,  and  eventually  calibration  must  be  relied  upon 
for  determining  this  response. 

An   estimate  of  this   behavior  may  be  obtained  from  the 
following  analysis  [Ref.  27], 

Consider  the  system  shown  in  Fig.  C-l.  The  fluctuating 
pressure  has  a  frequency,  ,  and  an  amplitude,  P  ,  and  is 
impressed  on  a  tube  of  length  L,  and  radius  r.  At  the  end  of 
this  tube  is  a  chamber  of  volume  V,  where  connection  with 
the  pressure  transducer  is  made  .  We  seek  an  expression  for 
the  pressure  indicated  by  the  transducer,  P  . 


V 


P0sinut 


T 


2r* 


% 


PRESSURE  TRANSDUCER 


Figure  C-l.   Schematic  Diagram  of  Pressure 
Transducer  Connection. 

The   mass  of  the  fluid   oscillates   under  the  influence 
of  fluid  friction  in  the  tube  which   tends  to  dampen  the  mo- 
tion. If  the  conventional  formula  for  laminar  friction  resi- 
stance in  tube  flow  is  used  to  represent  this  friction  ,  the 
resulting  expression  for  the  pressure  amplitude  ratio  is: 
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ft]  -  {[l  -<*>]»♦  ««,*>«}-! 


,-1/2 


(Eq.  C-l) 


where : 


GJ   is  the  natural  frequency  given  by  the  expression: 


2  2 
3  Trr  c 


%  " 


4LV 


(Eq.  C-2) 


and  h  is  the  damping  ratio  given  by  the  expression: 


h  = 


2Mg^  ,  Ulv 


(Eq.  C-3) 


per 


rr 


In  the  above  relations,  Eqs .  C-2  and  C-3  : 
c  =  the  velocity  of  sound  in  the  fluid. 
u  =  the  dynamic  viscosity  of  the  fluid, 
p  =  the  fluid  density. 
The  phase  angle  for  the  pressure  signal  is  given  by: 

.  ...-1  '  »"<-$;>  (B,.  0-4) 


+ 


'     (-ft;)1 


Introducing  the  physical  dimensions  of  the  present  appa- 
ratus, refering  to   Fig.  C-l   we  have: 

L   =    3,5    in  L   =    0.5    in 


r    =    0.25    in 


r' =    0.25    in 


or 


V    =nr2L    =    tt(  .25/12  )2  (3 .5/12 )    =    6.8177    10~5    (ft3) 

c    =    (K/p)1/2 

c  =  [144  (311  ,103)/1.938]1/2  =  4800  (ft/sec) 

p  =  62.34  lbm/ft3 

U  =  2.4  .10"5  lbf-sec/ft3   (at  60°F 
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Substituting  the  above  values  into  Eq.  C-2  we  get: 

3rr(. 25/12)2(4800)2 
UJn    4(3.  5/12M6.  8177  10"4) 


or:   u)   =  10885  sec 


-1 


Substituting  into   Eq.  C-3   we  get: 


2(2.4  10"5)(32.2)   |/3(3. 5/12)(6. 8177  10'4) 


h  = 


(62. 34)(4800)(. 25/12) 


TT 


-6 


or:    h  =  7.87  10 

From  these   results  we  see  that  the  difference   between 
pressure   transducer   indication   and  actual   plenum  chamber 
pressure  is  negligible  for  the  working  frequency  range. 

For  example  at  u)  =  1000  Hz  ,  from  Eq.  C-l   we  get: 


P 
P 


1- 


1000 


110885J    . 


+    4(7.87    10"6)2 


1000 


2-1/2 


10885 


or:         P/P      =    1.0085 
o 

Therefore  the  pressure  indication  error  due  to  dynamic 
response  is  0.85  % 

We  can  get  for  the  phase  angle  of  the  pressure  signal, 
from   Eq.  C-4  : 

4>    =    tan"1    -[2(7.87    lO-6 )(.  09187 )    /    (1    -    .09187)] 


or:       4>  =    -8.355    10 


-5 


rad 
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D.   CALIBRATION  DIAGRAMS 

In  order  to   compare  theoretical  and   experimental  re- 
sults concerning  the  indicated  pressure  from  the  mercury  ma- 
nometer, as  well  as  the  pressure  variations  from  the  digital 
multimeter,  it  was  necessary  to  calibrate  both  of  them. 
A  pressure  calibration  diagram  is  plotted  in  Fig.  D-l   where 
it  is  shown  to  be  linear. 

Since  the  Reynolds  number  is  a  function  of  flow  rate,  it 
was  easier  to  adjust  the  flow  more  accurately  with  the  aid 
of  a  pressure  indicating  digital  multimeter  than  the  rotame- 
ter. Therefore  two  diagrams  were  created,  shown  in  Figs.D-2 
and  D-3  where  are  plotted  the  indicated  flow  rate  versus  the 
indicated  pressure  ,  expressed  in  inches  Hg  and  D.C  Volts 
respect  i vely . 
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P(IN.HG)  VS  P(D.CVOLT) 
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Figure  D-l.   Pressure  Calibration  Diagram 
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FLOW  RATE  VS  P(IN.HG) 
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Figure   D-2.      Flow  Rate    vs.    Pressure    Indicated    (in.Hg)    Diagram 
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FLOW  RATE  VS  P(D.C  VOLT) 
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Figure  D-3.   Flow  rate  vs.  Pressure  Indicated  (V.DC)  Diagram 


94 


LIST  OF  REFERENCES 

1.  Rayleigh,  J.  W.  S.,  Philosophical  Magazine,  v.  34,  pp. 
59,  1892.  

2.  Tatsumi,  T.,"  Stability  of  the  Laminar  Inlet-flow  Prior 
to  the  Formation  of  Poiseuille  Regime, I, II  ,  "  Journal 
of  the  Physical  Society  of  Japan,  v.  7,  pp.  489-502, 
1952. 

3.  Batchelor,  G.  K.,M  Analysis  of  the  Stability  of  Ax i - 

symmetric  Jets  ,"  Journal  of  Fluid  Mechanics,  v.  14, 
pp.  529-551,  1962. 


4.    Crow,  S.  C. ,  Champagne,  F.  H.,  "  Orderly  Structure  in 
Jet   Turbulence  ,  "  Journal  of  Fluid  Mechanics,  v.  48, 
pp.  547-591,  1971. 


5.    Salwen,  H. ,  Grosch,  C.  E.,  "  The   Stability  of   Pois- 

seuille   Flow   in  a   Pipe  of  Circular  Cross-section  ," 
Journal  of  Fluid  Mechanics,  v.  54,  pp.  93-112,  1972. 


6.    Anderson,  R. ,  Bejan,  A.,"  The  Instability  on  a  Round 
Jet  Surrounded  by  an  Annular  Shear  Layer  ,  "  Journa 1 

of  Fluid  Engineering  ,  v.  107,  pp. 258-263,  1985. 


7.    Sarpkaya,  T. ,  "  A  Note  on  the  Stability  of  Developing 
Laminar  Pipe  Flow   Subjected  to  Axisymmetric  and  non 
Axisymmetric  Disturbances  ,  "  Journal  of  Fluid  Mecha- 
nics ,  v.  68,  pp.  345-351,  1975. 


8.  Garg,  V.  K.,  Rouleau,  W.  T. ,  "  Linear  Spatial  Stabi- 
lity of  Pipe  Poiseuille  Flow  ,"  Journal  of  Fluid 
Mechanics ,  v.  54,  pp.  113-127,  1972. 


9.    Leite,  R.  J.,"  An  Experimental   Investigation  of  the 
Stability  of  Poiseuille  Flow  ,"  Journal  of  Fluid  Mecha 
nics,  v.  5,  pp.  113-127,  1959. 


95 


10.  Darling,  G.  B.,  "  Heat  Transfer  to  Liquid  in  Intermit- 
tent Flow  ,  "  Petroleum  ,  v.  22  ,  pp.  177-180,  1959. 

11.  Sarpkaya,  T. ,  "  Experimental  Determination  of  the  Cri- 
tical Reynolds  Number  for  Pulsating   Poiseuille  Flow  ,  " 
Journal   of   Basic   Engineering,  (ASME),  v.  88,  pp. 
589-598,  1966. 

12.  Gilbrech,  D.  A.,  Combs,  G.  D.,  "  Critical  Reynolds  Num- 
bers for  Incompressible   Pulsating   Flow  in  a  Tube  , " 

Developments   in  Theoretical  and  Applied  Mechanics, 
Phenum  Press,  v.  1,  pp.  292-304,  1963. 

13.  Gilbrech  D.  A.,  Combs  G.  D.,  "  Pulsating   Flow  Re- 
search ,"  Uni  vers  i  ty  of  Arkansas . Eng ineer i ng  Experiment 
Station  Research  Paper  No  4  ,  1964. 


14.  Stuart,  J.  T. , "  On  the  Non-linear  Mechanics  of  Wave 
Disturbances  in  Stable  and  Unstable  Parallel  Flows  ," 
Journal  of  Fluid  Mechanics,  v.  9,  pp.  353-370,  1960. 


15.  Spielberg,  K. ,  Timan,  H. ,  "  On  Three  and   Two  Dimensio- 
nal Disturbances  of  Pipe  Flow  ,"  Journal   of   Applied 
Mechanics  ,  v.  3,  pp.  3  81-3  89,  1960. 

16.  Sparrow,  E.  M.,  Lin,  S.  H.,"  Flow  Development  in  the 
Hydrodynamic   Entrance   Region   of  Tubes   and   Ducts  ,  " 

The  Physics  of  Fluids,  v.  7,  pp.  338-347,  1964. 

17.  Klebanoff,  P.  S.,  Tidstrom,  K.  D.,  Sargent,  L.  M.  "  The 
Three  Dimensional  Nature  of   Boundary-layer  Instabili- 
ty ,  "  Journal  of  Fluid  Mechanics,  v.  12,  pp.  1-34,  1962 

18.  Miller,  J.  A.,  Fejer,  A.  A.,"  Transition  Phenomena   in 

Oscillating   Boundary-layer   Flows  ,"  Journal  of  Fluid 
Mechanics,  v.  18,  pp.  438-448,  1964. 

19.  Obremski,  H.  J.,  Fejer,  A.  A.,  "  Transition  in  Oscil- 
lating Boundary-layer  Flows  ,"  Journal  of  Fluid  Me- 
chanics ,  v.  29,  pp.  93-111,  1967. 

96 


20.  Landahl,  M.  T. ,  "  Wave  Mechanics  of   Breakdown  /'Journal 
of  Fluid  Mechanics ,  v.  56,  pp.  775-802,  1972. 

21.  Shen,  S.  F.,"  Some  Considerations  on  the  Laminar  Sta- 
bility  of   Time-dependant  Basic  Flows  ,  "  Journal  of 
Aerospace  Science,  v.  28,  pp.  397-417,  1961. 


22.  Grenspan,  H.  P . , M  On  Shear-Layer  instability  Breakdown 
and  Transition,  "  Journal  of  Fluid  Mechanics,  v.  15, 
pp.  135-153,  1963. 


23.   Uchida,  S.,  "  The   Pulsating   Viscous  Flow  Superposed 
on  the  Steady  Laminar   Motion   of  Incompressible   Fluid 
in   a  Circular   Pipe  ,  "  Zeitschrift  fur   Angewandte 
Mathematik  and  Mechanik  ,  v.  7,  pp.  403-422,  1956. 


24.   Hoyt,  J.  W.,  Taylor,  J.  J.,  "  Turbulence  Structure  in 
a  Water  Jet  Discharging  in  Air  , "  Physics  of  Fluids  , 
v.  20,  pp.  253-257,  1977. 


25.   Zucker,  R.  D.,  Fundamentals  of  Gas  Dynamics,  pp.  3A-1 
3D-6,  Matrix  Publishers,  1977. 


26.  White,  C.  M.,"  Streamline  Flow  Through  Curved  Pipes,  " 
Royal  Society  Proceedings  Series  A,  v.  123,  pp.  645- 
663,  1929. 


27.   Holman,  J.  P.,  Experimental  Methods  for  Engineers 
pp.  208-210,  Mc  Grow  Hill,  1984. 


97 


INITIAL  DISTRIBUTION  LIST 


No.  of  Copies 


1.  Defence  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia  22304-6145 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  California  93943 

3.  Department  Chairman,  Code  67  1 
Department  of  Aeronautics 

Naval  Postgraduate  School 
Monterey,  California  93943 

4.  Professor.  J. A.  Miller,  Code  67Mo  2 
Department  of  Aeronauutics 

Naval  Postgraduate  School 
Monterey,  California  93943 

5.  Hellenic  Air  Force  General  Staff  3 
ATTN:  Chief  of  Staff  HAF/GS,  Deputy  CS  HAF/GS 
Technical  Inspector  HAF/GS 

Holargos,  Athens,  GREECE 

6.  Hellenic  Air  Material  Command  1 
ATTN:  Com.  AMC 

Chief  of  Staff  AMC 
Elefsina,  GREECE 

7.  Air  Training  Command  2 
ATTN:  Com.  ATC 

Chief  of  Staff  ATC 
Tatoi-Dekel ia,  GREECE 

8.  Commander  Hellenic  Air  Force  Academy  5 
Tatoi-Dekel ia,  GREECE 

9.  Commander  Air  War  College  3 
Tatoi-Dekel ia,  GREECE 

10.  Hellenic  Air  Force  Research  &  Develoment  Center    4 
ATTN:  Com.  KETA 

KETA,  T.G.A  1010 

Pa  la  ion  Fa  1  iron,  Athens,  GREECE 

11.  Capt.  Z.Z.  Gikas  3 
4  Agorakritou  St.,  T.K  10446 

Athens,  GREECE 


98 


12.  Capt.  D.M.  Petridis 
SMC  1351 

Naval  Postgraduate  School 
Monterey,  California  93943 

13.  Lt.  E.S.  Mitrou 
SMC  27  7  7 

Naval  Postgraduate  School 
Monterey,  California  93943 


99 


/ 


211  2 


Gikas 

Effect  of  small 
pressure  disturbances 
on  the  breakdown  of 
round  laminar  and 
turbulent  jets. 


Thesis 
G4163 

c.l 


Gikas 

Effect  of  small 
pressure  disturbances 
on  the  breakdown  of 
round  laminar  and 
turbulent  jets. 


->e..,a> 


